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Diffusion's induced transport in periodic channels 

and an inverse problem 

Gershon Wolansk}0 



Abstract 



A diffusion's induced transport is defined for a linear model of a Fokker-Plank 
equation under periodic boundary conditions in one-dimensional geometry. The 
5J^. flow is generated by a diffusion and a periodic deriving force induced by a velocity 

C/^ I potential. An inverse problem is suggested for evaluating the deriving force in 

' terms of the response function associated with the flow. It is also shown that the 

inverse problem can be partially solved under some simplifying assumption. 

^' 

i '■ 1 Transport induced by a diffusion 

Recently, there is an increasing interest in transport process induced by an irreversible 
diffusion. This is motivated by some aspects of molecular dynamics and molecular 
motors [DEO]. 

^ ! The general paradigm of diffusion's induced transport is a model of several com- 

I ponents. Each component undergo an independent diffusion process, and the different 
^ I components are coupled together by the poisson process [SKB], [PJAP]. 

In some cases the basic model is that of a diffusion of a single component, while 

ON 
O 
00 

9: pt = {ap, + ^,p)^ a;G[0,l], (1.1) 



the transport mechanism is induced by an time periodic deriving potential [KK]: 



: (C^x + ^xp)iO,t) = (C^x + ^xp){l,t) = . (1.2) 

■ Here the deriving potential \E' := \E'(a;, t) = ^'(a;, t + T) is an appropriately defined 
function, and p > is the probability density of position of a "test particle". 

The model system (11. ip . as the one studied in [KK], deals with no-flux boundary 
conditions (11. 2p . It is proved that the solution of (ll.lfl.2p converges asymptotically, 
under general condition, to a periodic solution p{x,t) = p{x,t + T). The phenomenon 
of transport is introduced in a descriptive way: The period average of the asymptotic 
density p is shown to concentrate, under suitable conditions on \E', near one of the ends 
of the interval. 

The basic model introduced in [DEO], however, implies a geometry of periodic 
spatial structure. So, it seems reasonable to replace (11.21) by 

p(0,t) =p(l,t) , p,(0,t) =p,(l,t) . (1.3) 



^ Department of Mathematics, Technion, Haifa, Israel. E.mail gershonw@math.technion.ac.il 



1 



where \1/ is assumed to be 1— periodic in space (in addition to its time periodicity). 
Under these b.c, there is a natural definition for a diffusion transport, as follows: 

Consider a solution p = p{x,t) of (11.11 II. 3p . The velocity of a "particle" driven by 
this process is composed of the driving potential —"^x and the osmotic velocity —crpx/ p: 

v{x,t) = -^x- (Tpx/p (1.4) 

so the diffusion equation ( 11. ip takes the form of the continuity equation 

dp ^ djvp) ^ ^ 
dt dx 

Hence, we may describe the orbit x = x{t) of such a particle by 

^ = v{x{t),t) . (1.5) 
Recall that v is 1— periodic in x and T— periodic in t. The mean velocity, given by 

K := lim — ^ = lim - [ v{x{s), s)ds 

t^OO t t^QO t J 

is independent of its initial position x(0) and the initial data p(,0). It is just natural 
to define this process as a transport process if k 7^ 0. 

How can we prove that a certain driving potential \E' induces a transport k 7^ 0? 
As a case study we may consider a spatial periodic function %Ij{x) = tp{x + 1) and 

^{x,t) ■=ip{x-Vt) (1.6) 

for some \^ G M [applied voltage). This way we get a time periodic deriving potential 
of temporal period T = 1/V and spatial period 1. 

Switching now to the variable z = x — Vt, (11.11 11.31) takes the form 

Pt = {(ypz + '4'zP + Vp)^ , p{z, t) - p{z + 1, t) = p,{z, t) ~p,{z + l,t) = (1.7) 

where p{z,t) stands for p{x — Vt,t). In this case, the asymptotic limit p is just the 
stationary solution of (II. 7p which implies that the current / is a constant: 

ap, + ^,p + Vp = I . (1.8) 

This current is uniquely determined from V and the normalization condition 

lp=f Po = l. (1.9) 

From ([L5l OD and (Ol) 
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dz I 
dt p{z) 
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so from fll.9p 

z{t) 



lim 

t— >oc 

and 



V ~I 



Now, divide (11.81) by p and integrate over one period. Since both and ln(p) are 
periodic, we obtain 

K = i [\pr' 

Jo 

hence 



We now recall from (11.91) that /q^(p)^^ > 1, and J^iJ))^^ = 1 if and only if p = 1, which 
is the case if and only ii ip is a constant. We obtained 

Proposition 1.1. The process U.ll \1.3\) under condition U.6\) is a transport process 
( namely n 0) if and only if V ^ and ip is not a constant. Moreover, the mean 
velocity k is always in the direction of the applied voltage V, and \n\ < \V\. 



2 The inverse Problem 

In most applications we cannot measure directly the details of the driving potential ip. 
What we can measure is the current / as a function of the applied voltage V and the 
"temperature" a. 

The question we address now is related to the model problem (11.71) : 

What can be said about the driving potential ip from a complete knowledge of the mean 
velocity k(V, a) (equivalently, the response current /(y, cr)j, for each applied voltage 
V & and positive temperature cr > ? 

Indeed, it is shown in the appendix below that 

" \I{aV,a)J \ otherwise ^ ' ' 

where is the inverse Laplace transform action on functions of V into functions of 
z. 

So, a complete knowledge of the response current / as a function of cr and V yields 
a complete knowledge of the function 

F{a,z) := [ e'^'^^^'^-'^^'+'^^do; 
Jo 
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for any cr > and 5; G [0, 1]. I conjecture that a complete knowledge of F yields a 
complete knowledge on ip, up to and additive constant and shift. 

At this stage, however, I can only prove that a partial knowledge of /(V, a) yields, 
under certain assumption, a partial information on ip: 

At the limit = 0, the resistance ((9J/(9\^)^^(0, a) is given by 

Proposition 2.1. Assume ip is anti- symmetric ipi^x) = —ip^—x) (in addition to being 
1— periodic). Assume further we can measure the resistance {dI/dV)~^ as a function 
of the temperature a at zero voltage ^ = 0. Then we can find the distribution of ip. 

Indeed, if we know {dl /dV)y^^Q as a function of a > 0, then we can expand this 
function as a power series at a = cxd to obtain 

a-\dI/dV)~' = 1 + ^ + 4 + . . . . (2.3) 
Then, we expand the exponent on the left of (12.21) in powers of a^^: 

{ij{x) - ij{x + z)fdxdz (2.4) 



j e'^-'i^^-)-^^^+^))dxdz = l + J^^^ j j ( 



Recall that ip is 1— periodic. Comparing equal powers of a between (12. 3p and (12. 4p we 
get ci = and, for > 1, 



Cfc 

where Mjitp) = ijj^ the jth moment of ip. By anti-symmetry of ip we get Mj[ip) = 
for odd j. This implies, in particular, that Cfc = for odd k. In addition, the assumed 
known values of Ck for even k determine all even moments of ip, hence its distribution. 

Appendix 

We first calculate / = liV^a) from (II. 8p . Without limitation to generality we 
assume ip = Let 



a \ a 



and 



F±{z) := [ f±{s)ds . 
Jo 



The solution of fll.Sp takes the form 

where / and (3 are determined from the normahzation (11 .Qp and periodicity condition 
of p as follows: 

p(0) = p(l) =^ Pil - e-^/-) - Je-^/-F+(l) = (2.5) 

p = l ^pF.{l) + l[ UF+ = 1 (2.6) 
From (12. 5[ 12. 6p we factor out / to obtain 

^_ 1 - e-^/- 

" e-^/-F+(l)F_(l) + (1 - e-^/-) 

Integration by parts yields 



F+(1)F_(1)= f f-F^+ f UF_ . 
Jo Jo 



Let F+f_ = H+{V, a), = H_{V, a). Then 



Conversely, 



Let now 



e-''''^H.{V,a) + H^{V,a) = . (2.8) 

"-«^=1 • .^10.1] 

In particular (taking into account the periodicity of ip) 

h+{z + 1) + h.{z) = l[^i,o]{z) f e<^"'(^(^')-'^("+-))rfx (2.9) 

Jo 

where l[-i,o] is the indicator function of [—1,0]. We observe that 

/OO /'OO 
e-^''''h+{z)dz , H_{V) = / e''^'/"h_{z)dz. 
-OO J — OO 
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By (EZl ESD we get 



— oo 



oo 



e-^'/''{h4z) + h+{z + l))dz 



e^/- - 1 



So, if we take the inverse Laplace transform of (12.81) (as function of V/cr) we obtain 
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from ([21]) 




provided z G [—1,0], and C 




otherwise. 
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